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Abstract 

In this paper we consider the generation of a three-qubit GHZ-like thermal state by applying the 
entanglement swapping scheme of Zukowski et al. [Ann. N. Y. Acad. Sci. 755, 91 (1995)] to three 
pairs of two-qubit Heisenberg XY chains. The quality of the resulting three-qubit entanglement 
is studied by analyzing the teleportation fidelity, when it is used as a resource in the teleportation 
protocol of Karlsson et al. [Phys. Rev. A 58, 4394 (1998)]. We show that even though thermal noise 
in the original two-qubit states is amplified by the entanglement swapping process, we are still able 
to achieve nonclassical fidelities for the anisotropic Heisenberg XY chains at finitely higher and 
higher temperatures by adjusting the strengths of an external magnetic field. This has a positive 
implication on the solid-state realization of a quantum computer. 
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I. INTRODUCTION 

Algorithms such as quantum factoring Q and quantum search Q illustrate the great 
theoretical promise of quantum computers |3( . Consequently, many designs for the practical 
implementation of such devices have been proposed. Among these is one using quantum 
dots (localized electron spins) as quantum bits (qubits 4]), which has its roots in solid state 
physics j^ll^. The recent breakthroughs in the experimental physics of double quantum dot 
[3] show that this proposal is especially promising. Mathematically, a two-qubit Heisenberg 
XY chain at finite temperatures can serve as an effective model of a double quantum dot. 
The state of a two-qubit Heisenberg XY chain at thermal equilibrium (temperature T) is 
described by \ = e~P /Z, where H is the Hamiltonian, Z = tie~P H is the partition function, 
and 8 = 1/kT, where k is the Boltzmann's constant. The entanglement associated with the 

n 

thermal state y is referred to as the thermal entanglement I8IL 

Q 

In a very important 1999 paper, Ref. [9j, Gottesman and Chuang showed that single 
Qubit operations, Bell-basis measurements and Greenberger-Horne-Zeilinger (GHZ) states 

10j | are sufficient to construct a universal quantum computer. The generation of GHZ 
states therefore becomes an important problem. Slightly earlier, Zukowski et al. 11( and 
Bose et al. had shown how to obtain GHZ states from say "simpler" Bell pairs via 
entanglement swapping. Motivated by these experimental and theoretical developments, in 
this work, we consider as a first step the generation of a three-qubit GHZ-like thermal state 
from three pairs of two-qubit Heisenberg XY chains using the technique of entanglement 
swapping [lj]]. Since quantum teleportation [3| lies at the heart of Gottesman and Chuang's 
argument, we determine the quality of the resulting three-qubit entanglement by analyzing 
the teleportation fidelity, when it is used as a resource in the teleportation protocol of 
Xarlsson et al. [15], a generalization of the standard teleportation protocol of Bennett et al. 

"l\ to one using entangled three-qubit resource. 

This paper is thus organized as follows. In Section II, we give the Hamiltonian for a two- 
qubit Heisen berg X Y chain in the presence of an external magnetic field. Properties such as 



concurrence 



HQ and funded fr a Ctl0n Q Q of the — temal state a r e 



presented in Section III. In Sections IV and V, we briefly describe the entanglement swapping 
scheme of Zukowski et al. and the teleportation scheme of Karlsson et al. respectively. They 
are cast into forms suitable for our considerations here. We give the resulting three-qubit 
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thermal state in Section IV and the formula for calculating the correspoding teleportation 
fidelity in Section V. These set the stage necessary for the presentation of our results in 
Section VI. First, using the isotropic two-qubit Heisenberg XX chains, we establish that 
the entanglement swapping procedure amlpifies noise originally present in the two-qubit 
thermal states. Next, we show that despite this, one is still able to achieve nonclassical 
fidelities for the anisotropic Heisenberg XY chains at finitely higher and higher temperatures 
by adjusting the strengths of an external magnetic field. In the concluding Section VII, we 
summarize and discuss two possible research problems related to our results here. 



II. TWO QUBIT HEISENBERG XY CHAIN 

The Hamiltonian H for a two-qubit Heisenberg XY chain in an external magnetic field 
B m = 7] J (rj is a real number) along the z axis is 20] 



H = ~(1 + i)Ja\ ® a\ + ~(1 - ~i)Ja 2 A <g> a% + -B m {a\ ®a B + a A ® a 3 B ), (1) 

where <7° is the identity matrix and a l a (i = 1, 2, 3) are the Pauli matrices at site a = A, B. 
The parameter —1 < 7 < 1 measures the anisotropy of the system and equals for the 
isotropic XX modelj^lj and ±1 for the Ising model 2^|. (1 + 7) J and (1 — 7) J are real 
coupling constants for the spin interaction. The chain is said to be antiferromagnetic for 
J > and ferromagnetic for J < 0. 

The eigenvalues and eigenvectors of H are given by 

H\^) AB = £|$°> AB , 

H\&) AB = J\^)ab, 

H\$ 2 )ab = -J\$ 2 )ab, 

H\&) AB = -B\&) AB) (2) 



where B = J B^ + 7 2 J 2 = \A/ 2 + 7 2 J, 



= /,„ \, = [(B + B m )\00} AB + 1 J\ll) AB ] } (3) 
/ (i3 + 5 m ) 2 + 7 2 J 2 



\&)ab = -j=[\01)ab +\10) ab], (4) 
\&)ab = -±=[\01) AB -\10) AB ], (5) 



|$ 3 > AB = [{B-B m )\m) AB - 1 J\ll) AB \. (6) 

i{B-B m y + ^j* 



While l* 1 )^ = ^[\0l) ab + \W)ab] = I^bJab and \&) AB = ^[|01 /Afl ~ |10)ab] = 
\^Beii)AB are maximally entangled Bell states, |$°)ab and \§ 3 )ab are, so long as 7 ^ 
0, nonmaximally entangled states (see Eq.(18)). In fact, Eqs.(3) and (6) reduce to the 
maximally entangled Bell states -^[\00) AB + |11) AB ] = I^jWab and ^[|00) ab - |H)ab] = 
\^Beii)AB respectively, only when B m = (or 77 = 0). For 7 = 0, we have |$°)ab = ab 



and |$ 3 )ab = |H)abj product eigenstates with eigenvalues B m and — B m respectively [21 



III. THERMAL STATE, CONCURRENCE AND FULLY ENTANGLED FRAC- 
TION 

For the above system in thermal equilibrium at temperature T, its state is described by 
the density operator 

Xab = ^[e- /3B |$°)AB($°| + e- /3J |$ 1 ) AB ($ 1 | + e ^|$ 2 ) AB ($ 2 | + e^|$ 3 U($ 3 |], (7) 



where the partition function Z = 2cosh/3i3 + 2 cosh j3 J, the Boltzmann's constant k = 1 
from hereon, and (3 = 1/T. To quantify the amount of entanglement associated with xab, we 
consider the concurrence C[xab] = max{Ai — A2 — A3 — A4, 0}, where A& (k = 1, 2, 3, 4) 

are the square roots of the eigenvalues in decreasing order of magnitude of the spin-flipped 
density-matrix operator Rab = Xab{ca ® (7 b)Xab( 'a ® °b)> ^ ne asterisk indicates complex 
conjugation. After some straightforward algebra, we obtain 

Ai = \eP\ (8) 
A2 = -U"^ J , (9) 



A, 



Z 
1 



2^ 2 P 2^ 7 / 7 2 

1 + -^==- sinh 2 0B + -4- A/ 1 + -V" sinh 2 /3B sinh #B, (10) 



Aj - -m + ^^sinh 2 /3S--^-yi + ^sinh 2 /3Ssinh/3i3. (11) 

The concurrence derived from Eqs.(8) - (11) is invariant under the substitutions n — > —n, 
7 — ► —7, and J — ► —J. The same applies to the fully entangled fraction J-'Ixab]- 
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Therefore, we restrict our considerations to r\ > 0, < 7 < 1, and J > 24j : 



-^[xab] = max {ab^b^XabWbmI ab\ 

1=0,1,2,6 

\e^ J for v / j5 I +7 < 1, 

z , N ' ~ (12) 

I (cosh/?£ + $ smh/3#) for 07* + 7" > 1. 

For a bipartite entangled state pab to be useful for quantum teleportation we must have the 
fully entangled fraction F[pab] > |- Pab can then yield "nonclassical" teleportation fidelity 

$max (> 2/3) 

Eq.(7) reduces to the following three possibilities in the zero-temperature limit, i.e., 
(3 — > 00, at which the system is in its ground state. 



(a) if + 7 2 < 1: 

Xab = \[e pJ \&) AB {&\+eW\e) AB {e\] 

— > I$ 2 )^($ 2 |, (13) 

with Z = + e^. Eqs.(8) - (11) and (12) give C[xab] = 1 and F[x AB ] = 1 
respectively, in agreement with the fact that |$ 2 )ab is a maximally entangled Bell 
state. 



(b) 77 2 + 7 2 = l: 



2 

The above equally weighted mixture has 



XAB — i[|$ 2 )^($ 2 | + |$ 3 )AB($ 3 |]. (14) 



C[xab] = ^(1-7)t^0 (15) 



unless 7 = 1, but is useless for quantum teleportation since 



Hxab] = \ (16) 



(c) if + 7 2 > 1: 



Xab = \[e^\e) AB (^\+e^ J \^) AB (e\ 
,3 



&)ab(&\, (17) 
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where Z is as in Eq.(13). We have 



C[ X ab] = rz ^ z , (18) 



Vv 2 + i 2 

2 ( 1+ ^tT+tV " 2" 



^Ixab] = + ) > i ( 19 ) 



Therefore, 77 cr iticai = VI — 7 2 marks the point of quantum phase transition (phase transition 
taking place at zero temperature due to variation of interaction terms in the Hamiltonian 
of a system .§]) from the maximally entangled state Eq.(13), which yields C[xab] = 1 
and T\xab\ = 1, to the generally nonmaximally entangled state Eq.(17), which still yields 
nonzero C[xab] and !F{xab\ > | as long as 7 is nonzero and 77 finite: In the limit of large 77, 

C[ X ab] « IV' 1 — 0, (20) 

Hxab] « i + liv- 1 — i (21) 

when 77 is infinitely large. From Eqs.(15) and (18), we observe that for values of 7 other than 
7=|, there is a sudden increase or decrease in C[xab] at Critical, depending on whether 
7 > 3 or 7 < g, before decreasing to zero asymptotically, as 77 is increased beyond the critical 
value r/criticai- In contrast, -FfxAs] always increases with increasing 7 at //critical [see Eqs.(16) 
and (19)]. 

For nonzero temperatures, due to mixing of states, C[xab] decreases to zero and T[xab] 
decreases to | as the temperature T is increased beyond the critical values T^} ticaX and 
Critical respectively. Except when 7 = 0, both T^ tical and T^ tica i depend on 77. If 7 = 0, 
T^ icaJ = 1.13459J for all 77 [21], but T^ 2 ]^ decreases from 1.13459 J to zero as 77 is increased 
from to 77 cr iticai = 1, and remains zero when 77 is increased beyond 1 j^J. For nonzero 7, 
the behavior of T^ ticsl and T c ^ ical are qualitatively similar 2^. They both decrease with 77 



when 77 is increased from to 77 cr iticai = a/1 — 7 2 . However, in this case, this is followed by 



monotonic increase as 77 is increased beyond a/1 — 7 2 . Indeed, in the limit of large 77, 

T (2) V J (22) 

critical l nr? _ln 7 + l n2 - 1 ' 



Anisotropy therefore not only allows one to obtain entanglec 
and finitely higher B m than is possible in the isotropic case 



entanglement is useful as resource for teleportation 



24] 



cubits at finitely higher T 
J, but also the associated 
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IV. THE ENTANGLEMENT SWAPPING SCHEME OF ZUKOWSKI ET AL. 



1. (ill presented an explicit ent anglenient swapping scheme to generate a 



Zukowski et al. 

three-qubit GHZ state [10] from three Bell pairs:/^!?!, A 2 B 2 , A 3 B S . In their scheme, a local 
joint measurement in the three-qubit GHZ basis: 



I^GHZ) = 


V2 


|000) + 


1111)), 


\*GHZ) = 




|001> + 


|110>), 


\^Ihz) = 


V2 


|010) + 


|101>), 


\*GHZ) = 




|011) + 


|100», 


\*GHz) = 




|011) - 


|100», 


\*GHz) = 




|010) - 


|101», 


\*GHz) = 




|001) - 


|no», 


\^GHz) = 




|000) - 


1111)), 



(23) 

is performed on qubits Ai, A 2 and A3 from each of the three Bell pairs. They showed that 
this yields a GHZ state involving the other three qubits B\, B 2 and B 3 . Here, instead of 
Bell pairs, we consider xab, Eq.(7). The result is a three-qubit thermal state described by 



_ (*) 
XB 1 B 2 B 3 — Z^PiX-BiBiBa 



(24) 



8=1 



where 



Xb\b 2 B s = i-te AlA2A3 [(ll l AlA2A3 ® iB^B^iXAxBy ® Xa 2 b 2 ® Xa 3 b 3 ))- (25) 

Pi 



1 



The GHZ measurement is descibed by 

n^xAaAg = l^GHz) Ai A 2 A 3 (&GHz\ , ^ = 0, 1, 2, • ■ ■ , 7, 



and 



Pi = tr t( n A 1 A 2 A 3 ® lB 1 B 2 B a )(XA 1 B 1 ® XA 2 B 2 ® XA3B3) 



(26) 



(27) 



is the probability that the measurement outcome i occurs. The entanglement swapping 
scheme of Zukowski et al. was subsequently generalized by Bose et al. [12], to the case of 



7 



starting with N different multipartite GHZ states involving any number of qubits, doing 
local joint multipartite GHZ basis measurements by selecting any number of qubits from 
the N different multipartite GHZ states and also ending up with multipartite GHZ states 
involving any number of qubits. 



V. THE TELEPORTATION SCHEME OF KARLSSON ET AL. 



In Ref. |15( Karlsson et al. presented the teleportation of a two-level state \4>)a x using a 
three-qubit GHZ state, for instance \^ghz) a 2 bc, to either B or C in such a way that, only 
one of them, say C, can fully reconstruct the quantum state conditioned on the measurement 
outcomes of the joint Bell measurement on Ai and A 2 : 

n J Al A 2 = l*ktf>W*k«l, 3 = 0,1,2,3, (28) 
and the single-qubit measurement on B: 

U k B = \i; k ) B (i; k \, k = 1,2, (29) 

with = cos fi\0) b + sm fi\l) b and \ip 2 ) b = — sin//|0)£ + cos /x|1)b (0 < // < 71 /^)- Al- 

ternatively, the single-qubit measurement could be performed on C such that B reconstructs 
the quantum state. 

In this paper, reliability for teleportation will be the criterion for judging the quality of 
the resulting three-qubit state, Eq.(24). Therefore, in place of Y&ghz) a 2 bc we have x^a 2 bc 
[Eq.(25)]. The resulting state of qubit C after the measurements is then given by 

Pc 3k) = — tT AlA2 Bm Al A 2 ® n| ® i c )(|0)^(0| ® x%bc)1 (so) 

where 

q jk = tr[(ni lA2 ® H| ® 7c)(i0M0| ® xSUc)] (31) 
is the probability that the measurement outcomes j and /c occur. The teleportation protocol 
is completed by applying to p c an appropriate Pauli rotation, which depending on j and 
k, maximizes the teleportation fidelity: 

$max = fsinW X C(0|^ fc) |0)c, (32) 

where the integration is over all pure two-level states, |0) = cos 1 10) + e J</3 sin 1 1 1) (0 < i? < 7r, 
< </j < 27r). $ max quantitatively measures the reliability of Xa 2 bc as resource for the 
teleportation scheme of Karlsson et al. 
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VI. RESULTS 



For the two-qubit Heisenberg XY chains at finite temperatures [Eq.(7)], we have 

$max = Ci + C 2 cos fi sin n, (33) 

where 

_ 2(cosh 2 (3B + cosh (3B cosh (3 J + cosh 2 f3J) 
1_ 3(cosh/3£ + cosh/3J) 2 ' ^ ' 

_ 2 (sinh 3 (3 J + & sinh 2 f3 J sinh (313 + ^ sinh (3 J sinh 2 (3B + ^f- sinh 3 (3B) 
2 ~ 3(cosh/?S + cosh/5J) 3 P| (35) 

for both 77 2 + 7 2 < 1 and ?7 2 + 7 2 > 1. This is in contrast to C[xab] (see Ref. |20j]) and Eq.(12). 

In the zero temperature limit, Eq.(33) reduces to 

2 

$max = g(l + cos /i sin /i) (36) 



when ?7 2 + 7 2 < 1, a signature of (pure) three-qubit GHZ teleportation-resource (see Ref. 

11 2 

$max = - + — ( 1 + 7 + 7 2 + 7 3 )cos//sin//< - (37) 

when ?7 2 + 7 2 = 1, in agreement with Eq.(16); and 

\ 3 



*. 2 



max 



1+ WW+^) cos/iSin/i 



(38) 



when r/ 2 + 7 2 > 1. Therefore, $ max > 2/3 as long as 7 is nonzero and rj finite as in Eq.(21): 
In the limit of large 77, 

2 2 
$max ~ -(1 + 7V 3 cos/xsin/i) > -, (39) 

equals 2/3 asymptotically when 77 is infinitely large. 

For nonzero temperatures, we observe that both C\ and C2 are positive definite. Hence, 

for optimal $ max , we choose fj, — ir/4 from hereon. Due to mixing of states, $ max decreases 

to 2/3 at the critical temperature T cri ( ical beyond which the teleportation fidelity is worse 

(3) 

than what classical communication protocol can offer. For each 7, T cri ( ical is dependent on 
77. To obtain T cri £ ical we consider the following. For the thermal state Eq.(24) to be useful 
for quantum teleportation at nonzero T, we demand that 



7 



sinh d (3 J + ' = sinh 2 (3 J sinh pJr] 2 + 7 2 J 

v?r + 7 

7 V . , , , . , 2 , — - , ( 



[-^=1 sinh^Jsinh^V^ + ^+l-^^f] m^P^ + ^J 



> 2(cosh 2 (3 J cosh (3 J 3f + 7 2 J + cosh f3 J cosh 2 f3 J 3? + 7 2 J). (40) 



V 





0.1 


0.2 


0.3 


0.4 


0.5 


0.6 


0.7 


0.8 


0.9 


T (2) 1 J 

A critical / ° 


1.13459 


1.13105 


1.12029 


1.10193 


1.07525 


1.03904 


0.99126 


0.92828 


0.84267 


0.71411 


7^(3) / T 
^ critical/ J 


0.55508 


0.55093 


0.53825 


0.51631 


0.48371 


0.43810 


0.37605 


0.29473 


0.19890 


0.09950 



Table 1: Critical temperatures T c 5 ica j and T^ tical for various values of 77, when 7 = 0. 



When 7 = 0, Eq.(40) reduces to 



sinh 3 /? J > 2(cosh 2 /3Jcosh/??7J+ cosh/? J cosh 2 j3r)J). 



(41) 



Hence, T^ tical can be obtained from 



sinh 



T, 



"(cosh 2 — 77T, — cosh' 1 



(3) 



2 cosh" 



(3) 

critical 



,(3) 



cosh — 7^7 — cosh 2 ' 



t: 



(3) 



,(3) 



(42) 



■ critical \ critical critical critical critical ' 

It decreases from 0.555081J to as rj is increased from to Critical = 1, and remains zero 
when 7] is increased beyond 1 (see FIG. 1). This is because Eq.(41) is unattainable, consistent 
with the fact that the original two-qubit thermal states [Eq.(7)] have fully entangled fraction 
less or equal to 1/2 and thus cannot yield nonclassical teleportation fidelity, when 77 > 1. 
From Eq.(36) and a comparison 

(see Table 1), we conclude that the entanglement swapping process "amplifies" the thermal 
noise present in the original two-qubit states [Eq.(7)] resulting in T^ tical < ^critical ^ or eacn 
given 77. 

For nonzero 7, Eq.(40) can be satisfied as long as rf + 7 2 < 1 or rf + 7 2 > 1. At 
rf + 7 2 = 1, it reduces to 



1 + 7 + 7 2 + 7 3 ) sinh 3 (3 J > 4 cosh 3 (3 J, 



(43) 



which is again unattainable. Therefore, T^ ti j similarly decreases to zero when 77 is increased 



from to Critical = yl- 7^, as shown in FIG. 1. This is then followed by a monotonic 
increase in T^ tical as 77 is increased beyond VI — 7 2 - Indeed, in the limit of large 77: 



T (3) ~ 

critical ~ 



as in Eq.(22). The behavior of T c ^ ical is therefore qualitatively similar to T^ tical and T^ tical 



TjJ 



3 In 77 — 3 In 7 + In 2 ' 



(44) 



(i) 



,(2) 



(compare FIG. 4 in Ref. |20j and FIG. 1 in Ref. |24| with our FIG. 1). Hence, anisotropy 



allows one to obtain three-qubit GHZ-like thermal states Eq.(24), which yield nonclassical 
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FIG. 1: Temperature T^ tical / J at which the teleportation fidelity is less than 2/3, plotted as a 
function of 77 (magnetic field strength B m ) for various values of the anisotropy parameter: 7 = 
(solid line), 7 = 0.3 (dash-dotted line), 7 = 0.6 (thick dashed line), and 7 = 1 (thin dashed line). 
In each case, the teleportation fidelity is less than or equal to 2/3 in the region bounded by (and 
generally above) the relevant curve. Note that for any finite temperature, there is an r] for which 
the fidelity is strictly greater than 2/3. 

teleportation fidelities at finitely higher temperatures by applying an external magnetic field 
of finitely greater strengths. This is somewhat surprising in view of the following two facts. 
First, even in the zero temperature limit, \Q 3 )ab is far from "ideal": Eqs.(17) - (21). Second, 
from the above conclusion, one would naturally expect these imperfections to be amplified. 



VII. CONCLUSIONS 



In summary, we have shown that despite the fact our entanglement swapping procedure 
amplifies thermal noise, anisotropy in the original Heisenberg XY chains enables the result- 
ing three-qubit GHZ-like mixed state to yield nonclassical teleportation fidelities at finitely 
higher and higher temperatures by suitably increasing the strengths of an external magnetic 
field. In the light of the result of Gottesman and Chuang |9], our result has therefore a 
positive implication on the quantum-dot realization of a quantum computer. However, we 
should point out that although nonclassical, the teleportation fidelities are very close to 2/3, 
as can easily be deduced from Eq.(39). Hence, one immediate research problem would be 

n n 

to look at models like those considered by Anteneodo et al. [25] and Zhou et al. |26| to see 
if this situation could be improved. Next, multi-qubit GHZ states would be necessary in 
any realistic setting. To produce these states would require more "complicated" GHZ basis 
measurements. Therefore, another interesting research problem would be to study how the 
noise amplification scales with the complexity of these measurements. 
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